Cosmology in massive gravity 
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We argue that more cosmological solutions in massive gravity can be obtained if the metric tensor 
and the tensor E M „ defined by Stiickelberg fields take the homogeneous and isotropic form. The 
standard cosmology with matter and radiation dominations in the past can be recovered and ACDM 
model is easily obtained. The dynamical evolution of the universe is modified at very early times. 
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The discovery of accelerating expansion of the universe 
by the observations of Type la supernovae in 1998 [1, 
2] motivated the search for dark energy and modified 
gravity. Dvali, Gabadadze and Porrati proposed that 
general relativity is modified at the cosmological scale Q . 
It is also possible that graviton has a small mass which 
is comparable to the Hubble scale. The first attempt 
of a theory of gravity with massive graviton was made 
by Fierz and Pauli |4|. However, the linear theory with 
the Fierz-Pauli mass is in contradiction with solar system 
tests 0, Q . Recently, de Rham, Gabadadze and Tolley 
introduced a nonlinear theory of massive gravity Q that 
is free from Bouldware-Deser ghost [H, 0] . Cosmological 
solutions for massive gravity were sought in rJ- HI ■ The 
first homogenous and isotropic solution was found for 
spatially open universe in 11| and the massive graviton 
term is equivalent to a cosmological constant. The same 
solutions were then found for spatially open and closed 
universe in 12, [l3|. Langlois and Naruko then found 
more general cosmological solutions 16]. We follow the 
approach in [r| and proposed a new approach to get 
more general cosmological solutions. 

The theory of massive gravity is base on the following 
action Q 



Ml 
2 



/ d 4 xy^g~{R + mlU) + S„ 



(1) 



where m g is the mass of the graviton, the mass term 

U =U 2 + a 3 U 3 + a 4 U 4 , (2) 



The reference metric r\ a h is usually taken as the 
Minkowski one. The cosmological solution was first found 
for open universe and the mass term behaves like 
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an effective cosmological constant with 



A e ff = - m'g ( 1 + 3a 3 ± \j 1 + 3a 3 + 9a 2 , - 12a 4 ] x 



(1 + 9a§ - 24a 4 ± (1 + 3cv 3 )^/l + 3a 3 + 9aj - 12a 4 ) 



9(a 3 +4a 4 ) 2 



(8) 



The same solution was then found in 12[ for flat uni- 
verse. For a general case with positive, negative and zero 
curvature, Kobayashi et al. found the same solution with 
A e f f = m 2 ./ a for the particular choices of parameters a 3 
and a 4 [13|, 

a 3 = i(o!-l), a 4 = ^(a 2 - a + 1). (9) 

In [l6| , the authors assumed that the spacetime metric 
takes the form 



ds 2 = g^dx^dx" = ~N 2 {t)dt 2 + a 2 (t)-/ ij (x)dx i dx j , 



(10) 



with the spatial metric 



7-y {x)dx l dx 3 



dr 2 



1 — kr 2 



r 2 (dQ 2 



and generalized the reference metric from Minkowski 
metric to dc Sitter metric, 



U 2 = [K\ 2 -[K\ 
U 3 = [/C] 3 -3[/C][/C 2 ] +2[/C 3 ], 



and 



KS = SH- (VEX, 



(3) 
(4) 



Ui = [/C] 4 - 6[/C] 2 [/C 2 ] +8[/C 3 ][/C] - 6[/C 4 ], (5) 



(6) 



The tensor T,^ is defined by four Stiickelberg fields 4> a 
as 



Vabd^ a d^ b 



-dT 2 + bKT^jdX^X 3 



(11) 



where the Stiickelberg fields are assumed to be <p° = T = 
f(t), 4> l = X 1 = x l , so that the tensor E^^ takes the 
homogeneous and isotropic form, 



Diag{-/ 2 , bUfit)}^}, 



(12) 



and the functions bk{T) are 

bo(T) = e^ T , 6_ X (T) = H- 1 sinh(iT c T), 



(7) 



h(T) = H- 1 cosh(H c T), 
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they then found three branches of cosmological solutions, 
two of them correspond to the effective cosmological con- 
stant flS]) and exist for spatially flat, open and closed 
cases. They also found a new solution 

db k[f] a 

For the flat case, k = 0, substituting the de Sitter func- 
tion &o [/(*)] = e Hc ^ t%> into equation (fl"3|) . we obtain 
the effective energy density and pressure for the massive 
graviton, 

H 



Pg=~ m 2 g M% t \ l-^r 



H 2 

3(a 3 + ^ol 4 ) — 



H 



-3(1 + 5a 3 + 8a 4 ) — + 6 + 12a 3 + 12a 4 

tic 



(14) 



p g =m 2 n M 2 



3(a 3 + 4a 4 )— (l 



g pi 
6 + 12a 3 + 12a 4 



H 2 



(3 + 9a 3 + 12a 4 )-|- ( 3 + 



(15) 



H 2 



H 



(l + 6a 3 + 12a 4 ) — 3 + 2— 



So when H = H c , p g = 0. The Friedmann equations are 



H 2 + X = 



*M 2 pl 



(Prn + Pg), 



(16) 



2H + 3H 2 + - = ~— T (p m + Pg ). (17) 
a m pl 

For the flat case, substituting equations (|14|) and (fT5|) 
into Friedmann equations (|16|) and (|17l) . we get 



£ 2 (z) + 



H{Z) 



(a 3 + Aa 4 ) 



H 2 (z) 



(1 + 6a 3 + 12a 4 ) 



H{z) 



3(l + 3a 3 + 4a 4 ) 



(18) 



2-|(l + 2a 3 + 2a 4 ) 



ff 2 



3(1 + 3a 3 + 4a 4 )-^ - 2(1 + 6a 3 + l2a 4 )E(z) 
no 



+3(a 3 + 4a 4 )^E 2 {z) 
tic. 



-±E(z)-2E 2 (z) 



= 3n m (l + tv m )(l + z)^ 1+w -\ 

(19) 

where E(z) = H(z)/H$. The effective equation of state 
w g = Pg/ p g for the massive graviton is 

2E 2 {z)H/H 2 + 3Q m (l + w)(l + z)^ 1+w ^ 



w g = -1 



3[E 2 (z) - n m (l + z) 3 ^™)} 



Since p g = when -ff(z) = i? c , so if H c = ifn, we find 
that Q m — 1 which is inconsistent with current observa- 
tions, therefore H c ^ H . If H c < H , then we cannot 
recover the standard cosmology H 2 ~ p in the past un- 
less we fine tune the value of m 2 /H 2 to be very small. 
From equation (|T%1) , we see that the standard cosmology 
is recovered when H(z) <C H c and m g <C ifn < H c . At 
very early times, H(z) > H c , the universe evolves accord- 
ing to H 3 ~ p. If it was radiation dominated in the very 
early times, then the universe evolves faster as a(t) ~ t 3 / 4 
instead of t 1 / 2 . For the special case a 3 = a 4 = 0, Fried- 
mann equation is simplified to 




., . E 2 (z)-3 
EL„ I H c Hq 



3 —E{z)+2-% = Q m (l+z)^ 1+l 

tin 



At z = 0, -E(z) = 1, we get 



/7 2 



1 - Or, 



{H /H c -2)(H /H c -iy 



(21) 



(22) 



As discussed above, H c /H > 1, so to 2 , must be negative 
for this special case. Without loss of generality, we as- 
sume that to 2 = —j3iH 2 , and H c — with (3 2 > 1. 
For the special case a 3 = a 4 = 0, equation (j2"2"|) gives 



(1 - » m )/3 2 
(20,- 1)082 - 1)' 



(23) 



In this case, we have only two free parameters f2 TO and 
fa, and equation (flU)) gives 



_H_ 3» m (l + w m )(l + z) 3 ( 1+w "') 
H2 ~ 2(l-&l)E 2 (z) + 3%E(z) 



(24) 



If /?2 ^ 1, then /3i = (1 — Sl m ) /2, and the model becomes 
the ACDM model. This is shown in Fig. H]for /3 2 = 20.1 
and fl = 0.3. 

In summary, starting with the homogeneous and 
isotropic metric (fTt)|) and tensor (| 12[) , we obtain the Fried- 
mann equations (fT6| and §17\ with the effective energy 
density and pressure for the massive graviton, 



Pa 



Pg 



m 2 g M 2 



(b k [f}-a){6(l + 2a 3 + 2a 4 )a 2 



(3 + 15a 3 + 24a 4 )a6 fc [/] + 3(a 3 + 4a 4 )6 fc [/] 2 }, 

(25) 



to 2 M 2 i 



{ [6 + 12a 3 + 12a 4 - (3 + 9a 3 + 12a 4 )/]a 2 



2 [3 + 9a 3 + 12a 4 - (1 + 6a 3 + 12a 4 )f]ab k [f] 
[1 + 6a 3 + 12a 4 - 3(a 3 + 4a 4 )/]6 2 [/]}, 



(26) 



(20) and the equation of motion for the function f(t) which 
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FIG. 1: The evolution of the deceleration parameter q(z) and 
the effective equation of state w g of massive graviton. The 
blue lines are for the model with de Sitter metric as the refer- 
ence metric and the black lines are for the model taking &*,(/) 
as power law form. 



leads to the three branches of solutions 



bk[f(t)} 



(1 + 6a 3 + 12a 4 ± \/l + 3a 3 



9<*§ 



12a 4 ) 



3(ci3 + 4«4) 

db k [/] _ a_ 
df ~ N' 



a(t), 



(27) 
(28) 



When we take the solution (j2"7} . we get the effective cos- 
mological constant solution © independent of the choice 
of spatial curvature k. For k = 1, the solution was 
obtained in [ll| by taking the reference metric r/ a i, as 
Minkowski and the same homogeneous and isotropic ten- 
sor (H2]) with b k [f(t)} = f(t). The same solution ® 

I 



was obtained in 13j for all values of k for the particular 
parameters (|9|), but the tensor E MJ , is not homogeneous 
and isotropic. For the flat case k — 0, Gratia, Hu and 
Wyman obtained the same cosmological constant solu- 
tion by using another inhomogeneous and anisotropic 
tensor E MZ , 12j. However, the cosmological constant solu- 



tion (|8]) is just the consequence of the equation of motion 
of the function f(t) once we assumed the homogeneous 
and isotropic form for the metric (fT0| and the tensor T,^ 
(fl"2"|) . This suggests that the solution (|28|) is quite general 
even without the assumption of r/ a i, as de Sitter. There- 
fore, we propose that more solutions can be found from 
this approach by assuming more general form of &*[/(£)]. 

Follow the above argument, we assume a power law 
form b k [f(t)] = (-ff c /(i)) 7/(7_1) with 7 > 1, then the 
solution to equation (|28|) is 



&*[/(*)] 



0(7 - 1)H 



(29) 



and the effective energy density becomes 



Pg 



--3m 2 g M* 



3(1 



- (1 + 6a 3 + 
-2(l + 2a 3 - 

+(a 3 + 4a 4 ) 



3a 3 + 404) 
7-1 



12a 4 ) 
-2a 4 ) 
7- 



7 



II 



27-2 



7 



3 7 



37-3 



H 



,7-1 



H 



3 7 ' 



H 

K 

27 



(30) 



The effective pressure of massive graviton is 



Pg =m 2 Ml I 6(1 + 2a 3 + 2a 4 ) - 3(1 + 3a 3 + 4a 4 ) 



+ (1 + 6a 3 + 12a 4 ) 



7-1 



27 



,27-2 



-37(a 3 + 4a 4 ) 



7-1 



3 7 



a 



37-3 



H 



3 7 



27 



r 



l+2 7 1 



H 



H 2 



H 
H 2 



2 + 7 1 



H 
H 2 



(31) 



Again when the hubble parameter is in the range Hq < we require Ho <C H c . For the special case a 3 = 0:4 = 0, 
H (z) < i/ c , the standard cosmology is recovered. To Friedmann equations are 
have a long history of matter and radiation domination, 

I 



(32) 



4 



H 
H 



= -1 



Pi 



6/?l 



(33) 



with 



ft = 7 



(1 - ft m )^ 7 



(¥) 



2/3 2 7 



(34) 



In this case, we have three free parameters f2 m , 7 and 
/?2- Again if ,$2 ^> 1, pi = (1 — O m )/2 and the model is 
equivalent to the ACDM model. This is shown in Fig. Q] 
for (3 2 = 20.1 and fl = 0.3. 

In conclusion, more general cosmological solutions 
which are consistent with observational data can be 
found by taking homogeneous and isotropic form for both 
the metric and the tensor . The mass of graviton 
is in the order of ((1 - Q m )/2) 1 ^H Q . 
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